Evidence for coherent quantum phase-slips across a Josephson junction array 
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Superconducting order in a sufficiently narrow and infinitely long wire is destroyed at zero temper- 
ature by quantum fluctuations, which induce 2-n slips of the phase of the order parameter. However, 
in a finite-length wire coherent quantum phase-slips would manifest themselves simply as shifts of 
energy levels in the excitations spectrum of an electrical circuit incorporating this wire. The higher 
the phase-slips probability amplitude, the larger are the shifts. Phase-slips occurring at different 
locations along the wire interfere with each other. Due to the Aharonov-Casher effect, the resulting 
full amplitude of a phase-slip depends on the offset charges surrounding the wire. Slow temporal 
fluctuations of the offset charges make the phase-slips amplitudes random functions of time, and 
therefore turn energy levels shifts into linewidths. We experimentally observed this effect on a long 
Josephson junction array acting as a "slippery" wire. The slip-induced linewidths, despite being 
only of order 100 kHz, were resolved from the flux-dependent dephasing of the fluxonium qubit. 



I. INTRODUCTION 

The basic notion of superconductivity as a coopera- 
tive phenomenon is the superconducting order parame- 
ter. It is a continuous and complex single- valued function 
of coordinates, commonly referred to as the macroscopic 
wave function^. The electronic condensate of a super- 
conductor can move without friction, manifesting itself 
as a current without Joule dissipation. The density of 
a non-dissipative current (supercurrent) is proportional 
to the phase gradient of the macroscopic wave function; 
its amplitude is constant along the wire. At fixed gra- 
dient, phase is a linear function of the coordinate along 
the wire. For a long wire, the phase difference between 
the wire's ends may exceed, by many orders of magni- 
tude, the basic period 27r. This behavior is inconsistent 
with the thermodynamics of a superconductor, as its free 
energy is a 27r-periodic function of the phase difference. 
Thus, the equilibrium current reaches a maximum at a 
phase difference of tt and then oscillates with a further 
increase of phase difference. 

In a narrow wire (thinner than the coherence length), 
the adjustment of the supercurrent to the equilibrium 
value is achieved by the transient processes of 27r phase- 
slips (PS) of the macroscopic wave functiorP. As the 
wire becomes thinner, the phase-slips occur more fre- 
quently. Furthermore, at low temperature and for suffi- 
ciently small wire diameters, quantum phase-slips (QPS) 
take over the thermally activated ones, leading to an ac- 
tivationless relaxation of supercurrent^. 

It is important to notice that the conditions of the 
continuity and single- valuedness of the macroscopic wave 
function allow a 2nmi discontinuity of its phase at any 
point Xi along the wire (here are integers). At fixed 
phase difference between the ends of the wire, an observ- 
able quantity, such as current, depends on to = J2i '^ii 
but is independent of the specific locations Xi. This is 
why TO can be used to label different quantum states of 



the condensate. A ±27r-phase slip, occurring at some 
point along the wire, changes the value of to by ±1. Be- 
cause the state of the wire is characterized by to = J^i ™j 
rather than by each rrii separately, any of these ±27r-QPS 
results in the transition to — >■ to ± 1. The QPS processes 
a priori do not have to be dissipativ^^. In the absence 
of dissipation, QPS happening at different points along 
the wire interfere with each othei^I^. The resulting su- 
perpositions of QPS depend on the distribution of the 
electric charge along the wire due to Aharonov-Casher 
effeclP. We refer to these spatially interfereing QPS as 
coherent quantum phase-slips (CQPS). 

The thinner the wire, the larger the amplitude of 
CQPS and quantum uncertainty of to are. Proliferation 
of CQPS destroys superconductivity in ultra-thin long 
wires in the following sense: the equilibrium maximum 
supercurrent (prescribed by the periodic phase depen- 
dence of the free energy) decreases exponentially with 
the increasing length of the wir^, rather than being in- 
versely proportional to it. In general, CQPS are consid- 
ered to be the precursors of the superconductor-insulator 
quantum phase transitiorP. Moreover, in Josephson net- 
works with special symmetries, CQPS were predicted to 
give rise to the topologically non-trivial quantum collec- 
tive ordei'iSliil. On a practical side, challenging, but ulti- 
mately achievable control of external dissipation promises 
realization of a fundamental current standard by Bloch 
oscillations induced by CQPp^E^ES. 

In spite of their conceptual and practical importance, 
the CQPS processes have remained so far rather elu- 
sive. Several experiments have focused on the resistance 
of a current-biased superconducting wire, which there- 
fore dissipates an energy /<I>o per phase-slip event (J be- 
ing the bias current, $o = h/2e the flux quantum) ^ "^^ l 
Such experiments, while providing evidence for quantum- 
mechanical effects, cannot reveal the coherent aspects of 
quantum phase-slips and are likely to suffer from uncon- 
trolled dissipation in the electromagnetic environment 
of the biasing circuitry. Phase-biased chains of a few 
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Josephson junctions exhibit suppression of ma ximum su- 
percurrent consistent with the CQPS theorjff^l^, but 
have so far not displayed the expected spectrum of exci- 
tations associated with the quantum dynamics of phase- 
slips. Proposals for non-dissipative experiments with 
nanowires undergoing CQPS in the flux-qubit setufP^'^ni 
require a narrow range of quantum phase-slip amplitudes 
corresponding to transition frequencies of the order of a 
few GHz. Such experimental strategy appears feasible, 
but it may be complicated by the fact that phase-slip 
amplitudes of realistic wires vary by many orders of mag- 
nitude in view of their predicted exponential sensitivity 
to the wire transverse dimensions^. 

In this paper we describe an evidence for the coherent 
quantum phase-slips in a long array of A'^ = 43 Josephson 
junctions, which emulates a nanowire. Our experiment 
on CQPS features two important distinctions. First, for 
a faithful emulation of a nanowire, two parameters of the 
array are crucial: it must be long (comprising many junc- 
tions) and the probability amplitude of a phase-slip in a 
single junction must be small. These conditions minimize 
the effect of the discrete nature of the lumped circuit, 
and make the dynamics of condensate phase in the array 
similar to its counterpart in a wire. Second, to detect 
the presence of CQPS in our array, we take advantage of 
the sensitivity of the fluxonium qubit^ transitions to the 
phase-slips in its array inductance, which we treat as a 
slippery superconducting wire. 

This paper is organized as follows. In section|ll]of this 
paper we describe, both qualitatively and quantitatively, 
how the spectrum of fluxonium qubit is modified by the 
presence of CQPS in the qubit inductance. In section 
III we describe the experimental evidence of CQPS in 
the qubit inductance by comparing to the theory both 
frequency-domain and time-domain measurements of the 
fluxonium transitions. Appendix A provides details of 
our experimental techniques. Appendix B supports our 
theory and data analysis. 



II. DISPERSIVE EFFECT OF QUANTUM 
PHASE-SLIPS ON A SUPERCONDUCTING 
ARTIFICIAL ATOM 



and charging energy {Ecj = e'^/2Cj, Cj being junction 
capacitance). We neglect the capacitances of the array 
islands to ground (any large size metal object in the vicin- 
ity of the device). This is made possible by a tight pack- 
ing of array junctions (Fig. [I]i-inset), as justified in our 
previous wor]^^. 

For large Ej. /Eq^ , phase fluctuations are small. Cor- 
respondingly, quantum phase-slips crossing a ring made 
of such junctions occur rarely. This allows one to reduce 
the many-body Hamiltonian describing the quantum dy- 
namics of phases in the ring to a simplified, low-energy 
Hamiltonian^ . In the basis of states labeled by the num- 
ber m of phase-slips which crossed the ring the Hamilto- 
nian is given by (model A) 



H = ]^El{2^fY.^m - me.t)' H (m| 



E 



Eg \m) (m + 1| + h.c. 



(1) 



Here E]^ ~ (<I>o/27r)^ /Lx; corresponds to the inductive 
energy of the current in the loop of total inductance 
Ls = ^ Lj] the hybridization matrix element Eg is pro- 
portional to the CQPS probability amplitude in the loop, 
and we assume \E^\ <C Ef, so that one can neglect non- 
nearest-neighbor terms (like m o rn + 2); an externally 
applied magnetic field offsets the total flux in the loop 
by an amount <I>oxt = m-cxt'&o- 

In the absence of CQPS, i.e. Eg = 0, the energy lev- 
els of the ring depend quadratically on mext for a given 
m (Fig. [l];), and present two-fold degeneracy at half- 
integer TOcxt- The energy difference between the ground 
g and first excited e states varies in a zig-zag manner as a 
function of rricxt- CQPS removes the degeneracy, which 
leads to the rounding of the zig-zag corners at half-integer 
moxt- Moreover, the presence of CQPS allows e O g co- 
herent transitions under external radiation at any mext- 

Interestingly, the collective nature of the strongly- 
coupled junctions in this model hides in the specific form 
of the matrix element E^. Namely, 



N 



/2e 



(2) 



Fluxonium qubit (Fig. [T^, Fig. [5]) may be viewed as 
one junction shunted by the kinetic inductance of a se- 
ries array of A^ ^ 1 larger j unctions^ ^. Alternatively, we 
may view it as a loop of A^ + 1 junctions connected in 
series, with one "black-sheep" weaker junction (Fig. [T|d). 
The quantum state of the device is monitored by cou- 
pling the weak junction capacitively to an electromag- 
netic resonator , whi ch provides a non-dissipative mi- 
crowave readoutP2l2Sl. 

The loop array can first be understood qualitatively 
by the following semiclassical approach. Quantum fluc- 
tuations of the phase across a selected j-th Josephson 
junction are controlled by the ratio of its Josephson en- 
ergy {Ej. = ($o/27r)^/Lj, Lj being junction inductance) 



where Eg. is the "microscopic" contribution of an indi- 
vidual quantum phase-slip along the j-th junction'^ and 
Qj is the total charge on the islands between the 0-th and 

j 

j-th junction, i.e Qj — qi, with qj being the charge 

i=0 

on j-th island (Fig. ^p). As long as A^ ^ 1, one can 
view each quantum phase-slip event as tunneling of a fic- 
titious particle carrying flux $o from the inside of the 
loop to the outside (and vice- versa). Tunneling occurs 
via a superposition of multiple paths crossing different 
junctions and therefore encircle island charges. Tunnel- 
ing through each individual junction j then contributes to 
the total CQPS amplitude Eq. ^ with the correspond- 
ing Aharonov-Casher geometric phase^-^^^S^ 2TrQj/2e. 
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FIG. 1: (a) Optical image of the fluxonium device, bright color for Al, dark color for Si substrate. Two finger capacitors 
couple the small junction of the loop to the dispersive microwave reflectometry setup. Inset: Scanning electron micrograph of 
an array section, (b) Effective electrical circuit of the device with minimal elements for quantum phase-slip description. Each 
junction consists of a non-linear Josephson inductance Lj, a capacitance Cj, and an effective frequency-dependent resistance Rj 
representing possible intrinsic losses; every island carries offset charge qj induced by charged impurities in the oxide or unpaired 
quasiparticles. Pink arrows indicate possible semiclassical phase-slip trajectories. Charge-dependent energy shifts of artificial 
atom levels correspond to closed trajectories crossing several junctions (not represented on the figure), (c) Low-energy spectrum 
of the Josephson ring. Dashed parabolas represent energy of loop Josephson inductance threaded by an integer number m 
of flux quanta and tuned with external flux mcxt- (d) Two lowest transitions from the ground state, (e) Sketch of spatial 
distribution of superconducting phases of array islands, corresponding to the states with m = (no phase-slips) and m = 1 
(one phase-slip) across the entire loop at mext = 0. Phase of each island is represented by the angle between the corresponding 
arrow and the vertical. The value of m coinsides with the total phase winding along the loop devided by 27r. (f) Schematic 
broadening of the ground to first excited state transition <? — >■ e due to a spread of CQPS amplitude across the array. Origin 
of the spread is the combination of Aharonov-Casher effect and time-dependence of offset charges qj of array islands. Average 
transition frequency at mcxt = f/2 is dominated by the phase-slip across the weak junction. 



In real superconducting circuits, however, non- 
equilibrium charged impurities and quasiparticles cause 
the offset charges to fluctuate in time with an amplitude 
comparable to e, often on sub-ms time scale^^, much 
shorter than the averaging time of a typical qubit exper- 
iment. In the case of a homogeneous array (Eg. ~ Eg^ 
Vi,j), fluctuations of E^ and, as a result, fluctuations 
of the transition matrix element {g\m\e), which couples 



the qubit to the external radiation, are comparable to 
their respective means (see Appendix B 1| . This makes 
slow spectroscopic measurement at TOext = 1/2, and 
time-domain coherence measurements at any mext tech- 
nically challenging. We circumvent this obstacle by in- 
troducing an inhomogeneity in the array in the form of 
one {j = 0, for definiteness) weak "black-sheep" junc- 
tion. Since the phase across the weaker junction flue- 
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tuates more strongly, the corresponding phase-slip am- 
plitude Eso at j = largely exceeds that for all other 
(array) junctions, Esg ^ Esa = ^Sj^o^ where the 
averaging is taken over the junction index. Now, the 
Aharonov-Casher interference contrast is reduced from 
unity to a much smaller value of order N / 2Esa /Eso , 
where the \/N factor comes from averaging over ran- 
dom charges. Effectively, we split the roles of phase- 
slips in different junctions: while phase-slips across the 
black-sheep mix states of the loop with different m and 
therefore shape the qubit transition spectrum (Fig. ^p- 
c), the CQPS in the array junctions, combined with 
the fluctuating island charges, induce an inhomogeneous 
broadening Sv^g to the qubit g O e transition. This 
linewidth is maximal at |mext| = 1/2, where it is given 
by <5i^eg(l/2) = \/ N/2EsA/h, and diminishes away from 
that spot according to the following expression 



where 



^eg(l/2) 
l^eg(nT-ext) 



(3) 



(Fig. Ilf). The factor \/N/2 comes from averaging \E^\ 
over the random charges, assuming Esg ^ Esa- 

CQPS thus remarkably turn the flux-noise sweet spot 
into charge noise anti-sweet spot. We have now arrived 
at the main idea behind our experiment: by measuring 
the variation of the dephasing time T2 of the g ^ e 
transition of the fluxonium circuit as a function of ex- 
ternal flux moxt, we resolve the effect of CQPS in the 
Josephson junction array, provided that the transition 
intrinsic linewidth is smaller than the characteristic fre- 
quency y/N/2EsA/h of the CQPS. Evidence for CQPSs 
consists of: (i) the presence of the dephasing anti-sweet 
spot at |mext| = 1/2, (ii) excellent quantitative agree- 
ment of the observed dephasing time as a function of 
TOgxt with the described below parameter-free prediction 
Eq. Q , and (iii) confirmation of the inhomogeneous na- 
ture of the dephasing in echo measurements. The latter 
also allow us to extract limits for the time scale of charge 
re- arrangements. 

We introduce now a quantitative model of quantum 
phase-slips in the fluxonium circuit (model B). Unlike 
Eq. now we need to allow strong phase fluctuations 
in one of the junctions, while phase-slips in all others are 
still rare. Therefore, we use a mixed representation: the 
black-sheep junction is described by a continuous fluctu- 
ating phase ip; the dynamics of the remaining array of 
large junctions is treated by the same tight-binding type 
model as Eq. ([T]) in the space of the number m of CQPS 
across the array of large junctions only. The variables 
if and m are coupled because phase-slips in the array of 
large junctions create a phase bias on the black sheep. 
Thus, the Hamiltonian of model B is 



1 



E 



-Es \m) (m + 1| + h.c. 



(4) 



Hp (if, TOext) = -4£;c7^ - Ej cos ip 



+ -SL(27r)2((^/2^-mext)', (5) 

E,j = Ejg, and Ec = Ec^. The CQPS amplitude 
Es{Q} = E^{Q} - Es„ is different from that in Eq. ^ 
and accounts for phase-slips through every junction ex- 
cept the black-sheep one (we have set Qo = without the 
loss of generality). Similarly, the new inductive energy 
El excludes the contribution of the black-sheep junction 
inductance from the total inductive energy of the loop; 
I /El = I /El - l/Ejg. For Es = 0, expressions 
and ^ define the Hamiltonian of the fluxonium qubit!^. 
The second term in Eq Q incorporates the CQPS in 
the fluxonium inductance, and to now counts their num- 
ber. We are now in a position to evaluate accurately the 
flux dependence of the Aharonov-Casher linewidth of any 
a O /3 fluxonium transition. First order perturbation 
theory in \Es\ results in a remarkably simple expression 



(see Appendix B 1 for details of derivation) : 



EsA^/W,^ I 
OUap[m^^t) = |^a/3(mext)| 



(6) 



where 



Fapirric^t) ^ I d(^^'a(t/3)^'a((/3 - 27r) 
dc/?^'/3(<y9)*/3('y3- 27r) 



(7) 



with "ifaif) being the eigenfunction of the a-th energy 
state of the fluxonium Hamiltonian ([5]); it can be readily 
computed numerically. 

The dependence of Si^afi on the external flux mcxt 
comes from the presence of mcxt in Eq. ([5| determin- 
ing the wave functions ^a,^(v) which enter Eq. ([7|. For 
small amplitude of CQPS and for |moxt| close to 1/2, 
one recovers the relation ([3| from Eq. The overlap 
function Fap also appears in a previous worlP^ and can 
be understood from the fact that a phase-slip through 
the array inductance must shift the center of gravity 
of 5'q(i^) by 27r. The external flux dependence of the 
Aharonov-Casher linewidth thus encodes the overlap of 
the 27r-shifted fluxonium wave functions (Fig. [To| ) . The 
flux dependence of the linewidth (|6| of our model B gen- 
eralizes Eq. ([3]) of the intuitive model A to the case of 
arbitrary black-sheep junction parameters. 

Finally, since the finite lifetime of fluxonium excited 
states limit our experiment, and given the large number 
N of array junctions, a natural question arises: does the 
dissipation in the array scale up with A'^? Every junction 
of the array suffers from intrinsic, generally unknown dis- 
sipation sources, all lumped into a frequency-dependent 
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resistance Rj shunting the j-the junction (Fig.jljs). For- 
tunately (and counterintuitively) the dissipation in the 
black-sheep junction solely dominates the intrinsic re- 
laxation of the phase-slips spectrum with contribution 
from array junctions being suppressed as Indeed, 
the dominant phase-slip process generates a voltage pulse 
across the black-sheep; each array junction receives only 
a 1/iV portion of that voltage, and the total energy dissi- 
pated in the array resistors is only 1 /N of that dissipated 
in the black-sheep junction. We thus conclude that, re- 
markably, the apparent multi-junction complexity of the 
fluxonium circuit does not a priori penalize it with en- 
hanced energy relaxation. 



III. EXPERIMENTS 

In this section we present measurements of the transi- 
tion frequencies of a fluxonium qubit obtained using both 
frequency-domain and time-domain techniques, and an- 
alyze our data using the theory developed in section [ill 
We read out the qubit state with the help of a micro wave 
resonator, the technique known as circuit QEEEHHI. De- 
tails of this technique are described in the Appendix, 
Sections |A 2[ |A 3[ and |A 4| Capacitive coupling of the 
fluxonium circuit to the readout resonator results in the 
shift of the frequency of this resonator depending on the 
quantum state of fluxonium according to Eq. (A5). Es- 



sentially, the qubit state is mapped onto a particular 
value of the shift of the resonance frequency of the read- 
out resonator. This frequency shift is in turn detected by 



measuring the phase (Eq. (A2)) of the reflection ampli- 
tude (Eq. ( |Al| ) for a microwave signal, scattered off the 
resonator. 

Spectroscopy of fluxonium. Spectroscopy data in the 
range of qubit transition frequencies from 300 MHz to 
12 GHz and for the full span of external flux bias re- 
veals the spectrum of the transitions between the 3 low- 
est energy levels (Fig. [2]). The phase-slip frequency of 



the black-sheep junction. 



1/2) - \Es,\/h 



369 MHz corresponds to the center frequency of the line 
observed at Wext = 1/2 (Fig. [2^-left). By varying mcxt 
around that point and plotting spectroscopy traces on a 
2D color plot (Fig.[2]i-right) we observe the anti-crossing 
of the states of the loop with the phase-slips number m 
differing by a unity. Apart from the neighborhood of 
Woxt = 1/2 and m^xt = 0, the transition frequency de- 
pends linearly on the applied flux with a slope given by 
El, up to corrections of order El/Ej^ ~ 0.06 (Fig. [2]3). 
In niext — i^ plane (Fig.l2p), one recognizes the anticipated 
zig-zag shape (Fig. [ipj of the lowest g ^ e transition 

(Fig.d). 

The anticrossing at moxt = between transitions 
g f and 5 O e (Fig. [2j:) is associated with the hy- 
bridization between states m ^ m + 2. The size of this 
anticrossing according to model A would be orders of 
magnitude smaller than i^eg("^cxt = 1/2), the latter is 
associated with the hybridization m -H- m + 1. However, 



we find that the two frequencies are of the same order, 
because of the proximity of transitions 5 / and 5 O e 
at TOcxt = to the plasma resonance in the black-sheep 
junction, which model B properly takes into account. 

The 5 / transition anticrosses the cavity-assisted 
blue sideband of the g ^ e transition, which appears as a 
copy of the g ^ e transition (Fig.[2j: and inset). This an- 
ticrossing is an unusual vacuum Rabi resonance between 
the readout cavity and the e o / transition of the qubit, 
as indicated by the pole in the dispersive shift expression 
(A5) at v^f = vq. The magnitude of this anticrossing 
(Fig.[2j;-inset), which exceeds 100 MHz demonstrates the 
strong coupling of the qubit to the readout and accounts 
for the good qubit visibility in the remarkable 5-octave 
transition frequency range. The frequencies of all transi- 
tions are in perfect agreement with a numerical diagonal- 
ization of the hamiltonian of model B (neglecting Es). 
This allows us to extract accurately the circuit parame- 
ters values from spectroscopy data, see Table ( A 1 ). 

Measurement of decoherence times of fluxonium. We 
now turn to the time-domain experiments. The coher- 
ence of the g ^ e transition (Fig. [2k) is analyzed using 



standard time domain techniques in the full range of flux 
and frequency (—0.5 ^ Woxt ^ 0; 0.369 GHz ^ ^ 
9.1 GHz). First, the contribution of the energy relax- 
ation, is measured by the standard 7r-pulse techniques 
and yields the exponential decay times Ti of the excited 
state. Theoretical expression for Ti is given by Eqs. (B6 1, 
([BTI, pSl, and ^9^. In the narrow 300 MHz vicin- 



ity of the cavity frequency, Ti is limited by relaxation 
into the cavity (Purcell effect). Subtracting that effect, 
the overall shape of the Ti(TOext) data is matched well 
(within a factor of 2), by dissipation across the black- 
sheep junction through effective shunting resistance of 
the form R(){uj) = A/uj. Note that this agreement takes 
place for Ti varying by two order of magnitude and over 
a 4— octave frequency span. The only adjustable param- 
eter is the value oi A = 190 ± 60 MTl x GHz, the ex- 
treme values corresponding to top and bottom dashed 
lines (Fig. [3^) . Such 1 /a; frequency dependent resistance 
usually arises from the coupling to a large ensemble of 
discrete energy absorber^2ll^ whose microscopic origin is 
presently unclear'^. Dissipation across the black-sheep 
may likely come from dielectric losses in the coupling 
finger capacitors (Fig. [T^) . Alternatively, our energy re- 
laxation data could be explained by a lossless black-sheep 
and dissipation in the larger area junctions of the array, 
but with the A factor TV = 43 stronger. 

Dephasing times Tj* of the g ^ e transition, measured 
from the decay of Ramsey fringes, display pronounced 
minimum of about 250 ns at the flux sweet-spot and spec- 
tacularly rise by almost an order of magnitude (Fig. |3^- 
b), exceeding 2 /xs at |mcxt| — 0.2 {vg^ = 5.5 GHz). 
Around |moxt| = 1/2, the decay of Ramsey fringes is well 
fitted (Fig. [3}d) with a gaussian. This confirms the irrel- 
evance of the energy relaxation for Iwextl > 0.2, while it 
dominates at Woxt close to zero. Coherence times T2 ob- 
tained with echo experiments are drastically larger than 
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FIG. 2: (a) Left: Transition between g e fluxonium ground and first excited states at |moxt| = 1/2 measured by sweeping 
spectroscopy (Y-axis) frequency and measuring phase response of the 8.2 GHz microwave cavity coupled to the qubit (X-axis). 
Right: Raw spectroscopy data consisting of traces hke that on the left but with Y color-coded, as a function of external flux, 
in the vicinity of |mcxt| = 1/2; (b) Large scale spectroscopy data as a function of flux and frequency, (c) Extracted resonance 
locations vs spectroscopy frequency and applied flux. Solid line represent the fit to numerically diagonalized Hamiltonian 
neglecting quantum phase-slips in the array but taking into account interaction with readout cavity mode. Horizontal line 
represents cavity mode. Inset shows strong, 160 MHz anticrossing of the cavity-assisted blue sideband of the qubit lowest 
transition with the second lowest qubit transition, on which the readout mechanism is based. 



T2*, particularly around |moxt| = 1/2, and, for the most 
part, turn out to be limited by energy relaxation, i. e. 
T2 ~ 2Ti (Fig. [3^). Therefore, the noise that causes 
Ramsey fringes to decay is slow on the time scale of order 
T2 (about 10 lis) but fast on the time scale of the typi- 
cal Ramsey fringe acquisition time, of order one minute, 
typical of e-jump rates seen with sup ercon ducting single 
electron transistors and charge qubit^^ZIMJ, 

Experimental evidence for the CQPS. Analyzing our 
time-domain data, we find that the measured flux depen- 
dence of the decoherence times r2*("^ext) of the g ^ e 
transition is in excellent agreement with the expres- 
sions ^ and ([7|. Moreover, our data is inconsistent 
with conventional decoherence mechanisms, which in- 
clude low- frequency noise in the flux rriext , the Josephson 
energy Ej (due to critical current noise of the black-sheep 
junction), and in the inductive energy El (due to critical 
current noise in array junctions). 



tions in the array inductance, given by Eq. (B13), is also 



For the flux-noise mechanism, according to Eq. (Bll) 
and Fig. [2] we expect 2^2* ("^oxt) to be maximal at |r7ioxt| = 
1/2 (flux sweet-spot). This is obviously inconsistent with 
our data, which shows a clear minimum of T2*(mcxt) a-t 
I^Tiextl = 1/2, see Fig. [3^. Dephasing due to fiuctua- 



ruled out because the corresponding (mext) would also 
peak at |mext| = 1/2. (Qualitaively, this is because at 
1*710x11 = 1/2 the transition frequency v^g is given by 
\E^\/h, which is a property of the black-sheep junction, 
and is independent of the array inductance^.) The nu- 
merically evaluated sensitivity of the transition frequency 



'eg 



to El is presented in Fig. 



Theoretical analy- 



sis of the sensitivity of v^^g to the critical current noise 



in the black-sheep junction (Eq. ( |B12 1, Fig. 13 1) shows 
that the sensitivity must turn zero at some device-specific 
flux bias. Such non-monotonic dependence is inconsistent 
with the observed monotonic r|(mcxt)- 

Our key result is that the highly unusual flux- 
dependence of the dephasing time T2*(mcxt) (Fig. [3^) is 
well explained with the Aharonov-Casher effect of phase- 
slip interference. Upon adding the measured small con- 
tribution I/2T1 to the expression (B5) obtained from our 
theory of CQPS-induced dephasing, see Eqs. (|4|-([7]) and 
Appendix |B 1[ we find excellent agreement with the data 
(Fig. [3^). The only adjustable parameter we use is the 
theoretical value of at |TOoxt| = 1/2. Furthermore, 
this adjustable parameter agrees well with a WKB cal- 
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FIG. 3: (a) Summary of the coherence times measurements of the g e transition. Blue open circles denote energy decay 
times (Ti), red open triangles - decay time of the Ramsey fringes {T2) and closed green triangles the decay of the 7r-pulse echo 
experiment (T2). Black squares are the 7i data measured simultaneously with the echo data. The dash-dotted line represent 
the calculated relaxation time into the readout cavity (Purcell effect), which clearly is irrelevant for the most of the data. 
The two dotted lines are theoretical predictions of relaxation, for the two absolute values of effective resistance discussed in 
the text. Relaxation times Ti's present temporal variations that are not repeatable (differences between blue open and solid 
circles), while the overall dependence with applied flux is. Dashed line corresponds to the theory of the Aharonov-Casher effect 
induced dephasing, with the overlap function calculated numerically using circuit parameters extracted by spectroscopy, (b) 
Actual Ramsey fringes for several values of flux bias, with matched scale of the Y-axes for clarity, (c) Coherence time of the 
g f transition. Notations match those of (a), except that all relaxation times were multiplied by 2 for easy comparison with 
dephasing times. Theory lines obtained using same parameters as those for g o e transition. 



culation of Esa ^ hx 150 kHz based on our estimates of 
the array junction parameters. 

Control experiments. To verify our interpretation of 
the dephasing data in terms of the CQPS-induced broad- 
ening, we performed two main control experiments. 

In the first control experiment, we check if phase-slips 
interference can explain the dephasing of some other fiux- 
onium transition. For such experiment it is necessary to 
select a transition other than the lowest g ^ e, but with 
sufficiently long Ti , so that the decoherence is dominated 
by dephasing and not by energy relaxation. A good can- 
didate for this control experiment turns out the transition 
to the second excited state 5 o / in a narrow vicinity 
of rricxt = (Fig. [2j:). Interestingly, for this transition 
the flux bias mext — turns out to be a sweet-spot for 
Ti, see Fig. [Sj;. The sharp increase of Ti of the go/ 
transition at mext = is well explained, see Eq. ( |B10 ), 
by exactly the same model of dissipation, i?o(w) = A/ui, 
used to explain energy relaxation of the g ^ e transition, 
with A = 190 MQ x GHz. The origin of the peak in Ti 
at mext = comes from the fact that the / state can 
decay either directly to g or first to e and then to g; at 
"lext — 0, due to parity conservation, the direct f ^ g 
decay is forbiddcn^^. Thus, the bottleneck at mext — 
is a 280 MHz / •(-> e transition resulting in the lifetimes 
of the g f transition of up to 3 fj,s. 

The Ramsey fringe measurement on the g ^ f tran- 
sition at mext = 0, yields T2 ~ 2.5 /xs, a significantly 
smaller value than 2Ti ~ 6 /is. Echo measurement at 



TOext = yields a value of T2 almost matching 2Ti, 
see Fig. [3]:. Thus, there is a noticeable amount of in- 
homogeneous broadening of the 5 / transition at 
Woxt — 0. We find that this broadening is indeed pre- 
cisely accounted for by expressions Q, ([t]), and (B5|. 
In other words, the ratio of Tj* — 250 ns of the <? O e 
transition at |mext| = 1/2 and T2 = 2.5 /xs of the g ^ f 
transition at mext = coincides with the theoretical pre- 
diction Fgf{0)/Fge{l/2), where the overlap integrals are 
computed without any adjustable parameters. The re- 
duction of T2 with the departure from the sweet-spot 
Wcxt = is also correctly reproduced without adjustable 
parameters (Fig. [sj:). 

The second control experiment is even more telling, 
and is performed on another fluxonium device, accord- 
ing to the following logic. Our theory of dephasing by 
the interfereing quantum phase-slips predicts that in- 
creasing the Josephson energy Ej.^^ of the array junc- 
tions while leaving their charging energy Eq^^o the same, 
would suppress exponentiall}'^ the phase-slip matrix el- 
ement Esa- (Let us remind that we assume all the junc- 
tions of the fluxonium array to be nearly identical, i.e. 
Ej. ~ Ej^ and Eq^ — Ej. for any i, j > 0, so that 
\Esi\ — \Esj \ — Esa-) Thus, we shall be able to prac- 
tically switch off the decoherence of fluxoium transitions 
due to CQPS in its array inductance, by a small adjust- 
ment of the array junctions parameters. 

Parameters of the fluxonium device for the second 
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FIG. 4: Measurement of the dephasing times T2 of a new flux- 
onium device with supressed CQPS rate (magenta inverted 
triangles) and theoretical prediction (magenta dotted line). 
Data from Fig. 3a of the main text (red triangles) on T2 of the 
main device is shown for easy comparison. For |mcxt| < 0.25, 
energy relaxation dominates the decoherence. 



control experiment were extracted from spectroscopy 
data (not shown) to be Ej^^^ = Ej = 12.0 GHz, 
Ec^^o = Ec = 2.46 GHz, El = 0.89 GHz, very simi- 
lar to the parameters of the main device (see Table All. 
Since we kept the geometry of the array junctions and 
their number identical to the previous device, we infer 
that the charging energy of the array junctions Ec^^„ 
did not change significantly, while the Josephson energy 
Ej^^„ ^ N X El ^ 43 X 0.89 = 38 GHz increased by 
a factor of 1.7. For such device parameters, the spec- 
trum of this fluxonium is similar to that of the main one. 
However, in sharp contrast to the main device, we find 
that T2*(moxt) for the g -H- e transition is now sharply 
peaked at |mcxt| = 1/2 and is nearly independent on 
TOoxt away from this spot (Fig. |4]). The reduction of 
T2 away from the sweet-spot follows the prediction of 



the first-order flux-noise effect, Eq. (Bll), assuming the 
flux noise amplitude 5moxt ~ 10"'^. The dependence 
T2*(TOext) away from the spot |mcxt| = 1/2 is now clearly 
inconsistent with the CQPS. (They still may be respon- 
sible for the dephasing at |moxt| = 1/2.) The measured 
value of r| =4 ^s exactly at |mext| = 1/2 is too small to 
be explained by the second order flux-noise, but matches, 
within the factor of two, the prediction of Eq. ^ for 
the CQPS-induced dephasing. Thus, a second fluxonium 
device with nearly similar parameters, apart from an in- 
crease in Ej.^g by a factor of 1.7 resulted in a 16-fold 
enhancement of the coherence time at |TOoxt| = 1/2, in- 
dicating the expectecP^ exponential suppression of the 
phase-slip interference. 



IV. SUMMARY 



Coherent quantum phase-slips hybridize states of a 
superconductor with different configurations of the or- 
der parameter. The corresponding hybridization energy 
is exponentially small for "good" superconductors. In 
our experiment this energy, \Es\, corresponds to only a 
few hundred kHz, three orders of magnitude lower than 
the 15 mK sample temperature, and five orders of mag- 
nitude lower than the main qubit energy scales, i. g. 
plasma frequency and the inductive energy (both of or- 
der 10 GHz). Detection of the tiny energy scale \Es\ 
has been made possible by two ingredients specific to our 
experiment. First, the Aharonov-Casher modulation of 
CQPS amplitude broadens the qubit transitions. Second, 
the the immunity of the fluxonium circuit to dissipation 
(Q = I'egTi > 10^) allows us the high-precision measure- 
ment of this broadening, thus revealing CQPS. By replac- 
ing the junction array of the present experiment with an 
amorphous superconducting nanowire, but keeping the 
black-sheep junction, one may attempt seeing CQPS in 
nominally continuous wires with poorly controlled CQPS 
amplitudes. 

Our experiment also shows that the fluxonium artifl- 
cial atom may flnd applications in various quantum infor- 
mation processing schemes: it provides a 3-level system 
displaying a combination of larger frequency range and 
anharmonicity than most other qubits; it can operate 
away from flux sweet spots without loosing too much co- 
herence; its coupling to a microwave cavity can be varied 
from weak to strong for exchange of quantum informa- 
tion without the side effect of spontaneous emission. In 
the event that critical current noise would end up domi- 
nating superconducting qubit coherence, one may expect 
a 1/y/N suppression of this effect using an N junction 
array. Finally, demonstrated coherence quality factor 
Q exceeding 10^ in a circuit containing as many as 44 
junctions encourages the design of large quantum Joseph- 
son networks, especially those offering topologically pro- 
tected ground states. 
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Appendix A: Experimental techniques 
1. Sample description 

The qubit and the readout circuits (Fig. [5^) are fab- 
ricated on a Si chip using Al double-angle evaporation 
through a suspended electron beam resist mask. The 
readout part consists of a resonator, implemented with a 
A/4 coplanar-strips (CPS) transmission line. Resonator 
is coupled to 50 n measurement leads (microstrips) us- 
ing a pair of interdigitated capacitors (Fig. ^) . The two 
leads of the "black-sheep" junction are connected to the 
resonator strips with another pair of smaller interdigi- 
tated capacitors (Fig. [sj;) . A number of resistance and 
dose test structures are placed outside the resonator. 

The fabrication procedure is outlined in the previous 
worlPi], with relevant device parameters collected into 
Table (All. We emphasize that the entire fluxonium 



device - including the Josephson array, the black-sheep 
junction, the resonator, and the test structures - is fabri- 
cated in a single step of e-beam lithography and double- 
angle evaporation. Such simplification in the fabrication 
of a superconducting qubit has been made possible be- 
cause: i) the dimensions of both the black-sheep junction 
and the array junctions are chosen to be sufhciently close 
for patterning both types of junction in a single resist 
mask ii) the strips of the CPS resonator are sufficiently 
narrow so that the e-beam and lift off process used for 
the small junction fabrication could be readily applied to 
the resonator; in addition, space remains to accommo- 
date test junctions and arrays fabricated simultaneously 
with the qubit. 



2. Microwave reflectometry readout 

We use the lowest differential mode of our resonator, 
which corresponds to the frequency ujq at which the phys- 
ical length of the CPS transmission line matches a quar- 
ter of the wavelength. The next order resonance lies at 
3luq. The qubit can be viewed as a high impedance ter- 
mination Zl">(w) placed at the open end of the trans- 
mission line, and depends both on the frequency uj and 
the qubit state a (Fig. 6^) . It provides a small contribu- 
tion Xa ~ termed the dispersive shift - to the resonance 
frequency ujq, when the qubit is in state \a). The shift 
in the resonance frequency is detected by monitoring the 
complex-valued reflection amplitude P for the scattering 
of a microwave signal off the resonator. Approximating 
the single-mode resonance with an effective LC-oscillator 
(Fig. |6]d), the reflection amplitude is given by 



r(w,wo) 



2i{'- 



Qcxt + Qi, 



2^{^) + Q:A + Q^r 



(Al) 



where Qoxt is the quality factor due to the energy loss in 
the matched 50 measurement leads while Qint is the 
quality factor due to the energy loss inside the resonator 



(electrically represented by a resistor shunting the LC- 
circuit). In our CPS resonator Qoxt — 400 <C Qint ~ 
4000 making |r(a;, wg)! to be very close to unity (Fig.[6j;). 
The phase 6 — Arg(P) of the reflected signal is a rapid 



function of frequency: 9 = 2 arctan 2Qc 



(Fig. 6h 



Finally, the difference in phase of the reflection coefficient 
^eg = de — 9g between the qubit excited state e and the 
ground state g is given by 



= 2 arctan 2Q, 



2 arctan 2Q, 



(A2) 

where = cjo/Stt is the resonator frequency in Hz. The 
quantity 9eg is plotted on the Y-axes of Fig.[2]i-left, and 
plotted as the color scale in Fig. [2^-right and in Fig. [2]:. 

Note that since 9eg saturates quickly once 
2Qcxt ^"^^^ > 1; increasing Qcxt above the value 
wo/IXe ~ Xg\ '^ill not improve the sensitivity. In our 
experiment } ~ Xg\ li^ between 1 MHz and 10 MHz, 
while vq = 8.175 GHz, so Qcxt = 400 is a convenient 
choice with a measurement bandwidth of 20 MHz. We 
discuss the origin of the dispersive shifts as well as the 
effect of finite Qoxt on the qubit lifetimes (Purcell effect) 
in the text below, see Eqs. (B6), (B7|, and (B8). 



3. Measurement schematic 

Low-temperature setup (Fig. [7^). The incoming sig- 
nals line is attenuated with cryogenic high-frequency re- 
sistive film attenuators (XMA), with total attenuation 
exceeding 50 dB. The readout line is shielded by the 
two 4—12 GHz isolators (Pamtech) and a low-pass fil- 
ter (K&L) with rejection band 10 — 40 GHz. Outgoing 
signals are amplified using a 5 K noise temperature cryo- 
genic HEMT amplifier (Caltech). Incoming and outgoing 
waves are separated from each other using a directional 
coupler (Krytar). Differential excitation of the resonator 
is implemented using a 180 degree hybrid coupler (Kry- 
tar) . The chip rests at the bottom of the fully enclosing 
copper sample holder, the microwaves are guided to it 
by means of two printed circuit board microstrips wire- 
bonded to their on-chip continuation, and perpendicular 
coaxial-to-microstrip transition implemented using An- 
ritsu K connectors (Fig. [8|. The transition provides less 
than 15 dB return loss up to 20 GHz. Flux bias is 
provided by driving a hand-made superconducting coil 
(glued to the sample holder) with DC currents of order 
1-10 mA. 

Room-temperature setup (Fig. (tJd). The readout sig- 
nal is provided with Agilent E8257D generator (RF), the 
qubit pulses are generated using Agilent E8267D vector 
signal generator (IQ) combined with Tektronix 520 ar- 
bitrary waveform generator (AWG). Both readout and 
qubit signals are combined at room temperature and 
sent into the IN line of the refrigerator. The reflected 
~ 8 GHz readout signal from the refrigerator OUT line 
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FIG. 5; (a) Optical image of the sample (color filters applied for better contrast), bright indicates Al, dark indicates Si substrate. 
Sample contains the CPS resonator, the qubit, measurement leads and test structures, (b) Zoom in on the voltage antinode 
region of the resonator, (c) Further zoom on the fluxonium qubit loop. 
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FIG. 6: (a) Electrical model of a qubit interacting with the CPS transmission line resonator. Thick solid line represents the 
distributed transmission line, (b) Simplified circuit model of the dispersive effect of qubit on resonator, (c) Modulus of the 
reflection amplitude in model of panel (b) vs. signal frequency for Q^^t = lOQcxt- Blue and red traces correspond to qubit 
states g and e respectively, (d) Phase of the reflection amplitude, same conditions as in panel (c). 



is amplified at room temperature with a Miteq amplifier 
(1 — 12 GHz, 30 dB gain), mixed down with a local oscil- 
lator signal (LO), provided by HP 8672A, to a 50 MHz 
IF signal, then filtered and amplified with the IF ampli- 
fier (SRS SR445A), and finally digitized using one chan- 
nel of the 1 GS/s Agilent Acqiris digitizer. A reference 



IF signal is created by mixing a copy of RF and LO and 
digitized using the second channel. A software procedure 
then subtracts the phases of the two IF signals, resulting 
in a good long-term stability of the phase measurement. 
The short-term stability is implemented by phase lock- 
ing every instrument to a Rb 10 MHz reference (SRS 
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FIG. 7: (a) Dilution refridgerator part of the microwave reflectometry setup. HEMT - high electron mobility cryogenic amplifier, 
ISO - cryogenic broadband isolator, LP - low pass filter, DRC - directional coupler, HBC - 180 degree hybrid coupler, see text 
for more details (b) Room-temperature signal generation and demodulation setup. LO - continuous microwave source, RF 
- pulsed (Mk.l controlled) wave source, IQ is a vector signal generator for driving the qubit transitions, AWG - arbitrary 
waveform generator, A2D- two-channel fast digitizer, - matched signal combiner/splitter, "x" - microwave mixer, BP- and 
pass filter, "10 MHz" - Rb reference, RF AMP - a room temperature microwave amplifier. 




FIG. 8: Photographs of the chip holder carrying the actual fluxonium device. Anritsu K connector guides microwaves from 
coaxial cables to Cu microstrip lines (300 fim wide) on the printed circuit board. The center conductor of the K connector is 
soldered to the begining of the microstrip. Microstrip lines continue to the chip by means of 5-6 short wirebonds. 
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FS725). The marker signals of the AWG are used as 
triggers to other instruments. Typical time to acquire 
Ramsey fringes of 5 fjs long (Fig.jsjj) ranges from 10 sec- 
onds to 1 minute, without noticeable change in the fringe 
decay time. The magnetic coil is biased with Yokogawa 
7751 voltage source in series with a 1 : 10 voltage divider 
and a 1 kfi resistor at room temperature. 



Dispersive shifts of a resonator by fluxonium 
artificial atom 



to the / state involve large charge motion and the Vef 
frequency remains in a window 1 — 2 GHz away from the 
cavity frequency vq for all flux values. This behavior is 
due to a large participation of the black-sheep plasma 
mode in the / state for the present device parameters. 
By contrast, the g O e transition involves little charge 
motion (Fig. [9]) because it connects states with different 
phase-slips number m, and, in addition, detunes quickly 
from the cavity with the external flux. The large splitting 
between the transition e o / and the resonator shown in 
the inset of Fig.[2j; of the main text illustrates this point. 



Here we evaluate Xa ^^sed in Eq. (A2| for the reflec- 
tometry signal (Fig. [6|l). The interaction of the qubit 
with the resonator is given by 



V = hg{a + d^)n , 



(A3) 



where n is the charge on black-sheep junction capaci- 
tor in units of 2e, and a is the annihilation operator for 
the microwave photons in the equivalent LC-oscillator 
(Fig. [6^-b). The coupling constant g, given by 



Cj + Cc 



-Zq/Rqi'o 



(A4) 



is expressed via the characteristic impedance of the os- 
cillator Zq = AZoo/t^ {Zrxi — 70 — 80 O is the wave 
impedance of the CPS transmission line), and the super- 
conducting impedance quantum Rq = h/{2e)^ ~ 1 kO. 
Treating y as a perturbation to second order in g yields 
the following expression for the dispersive shift Xa the 
resonator frequency for the qubit in state \a): 



,,2 



(A5) 



where n^^ are the matrix element of the charge oper- 
ator fi connecting states a and j3 and is the qubit 
transition frequency (in Hz) between these states. In for- 
mula Vaj3 is taken negative if the state a is higher in 
energy than the state /3, and positive otherwise; in the 
rest of the text, for simplicity, we treat as a positive 
number. Perturbation theory breaks down whenever the 
denominator goes to zero, a situation which corresponds 
to a resonance between various qubit transitions and the 
resonator frequency. The inset of Fig. [2]; shows an un- 
common instance of such vacuum Rabi resonance with 
the qubit transition e o /. A conventional vacuum Rabi 
resonance, involving the lowest (7 o e qubit transition. 



takes place at |moxt| — 0.05 and 



.2 GHz, it is 



shown as an anticrossing of the theory lines (actual data 
available elsewhere^). 

Interestingly, for the case of the present fluxonium arti- 
ficial atom, the dominant contribution to the dispersive 
shift of the lowest transition {g o e) comes from the 
transitions g ^ f and e o / and not from the g e 
transition itself. This happens because the transitions 



Appendix B: Decoherence of fluxonium 

1. Aharonov-Casher Line Broadening 

Here we derive and discuss Eq. (l6|. The Hamiltonian 
of model B, defined by Eqs. Q and ([5| is invariant under 
the transformation (tp, m) — >■ (ip — 2Tr, rh+1). This sym- 
metry represents the fact that by looking only at the ini- 
tial and final states of the junction loop we cannot distin- 
guish which part of the loop (black-sheep junction or the 
array) actually underwent a phase-slip. This point can 
also be illustrated by representing graphically the phase 
distribution of array islands before and after a phase-slip 
by 27r(Fig.[l]3). The unperturbed eigenstates | a) of the 
Hamiltonian ^ then take the following form: 



\(o) _ 



lim , 



M 

^ E 

rn=-M 



^a(<y5 ~ 27r?7i) |m) 



(Bl) 

Here ^a('p) is the wavefunction of the a-th (non- 
degenerate) eigenstate of the fluxonium Hamiltonian 
Eq. ([5]), the states |m) are the eigenstates of the integer 
fa operator, and the normalization is chosen to satisfy 
(a|/3)*-°'' = 5ai3- Now, treating the quantum phase-slips 
perturbation (second term of Eq. Q) to the first order in 
I-EsIj we find that the correction to the qubit transition 
frequency v^^^ between the states a and 13: 



lim 



2h M',M"^oo ^2M' + W2M" + 1 



+ 00 



M' 



M" 



171— — 00 ni' — — M 



E E 



Alp 



X {^a{v - 2T:m')^a{v - 27rm") 
- ^[j{Lp - 2mji')^i3{(p - 2Trm")) 
X {Eg (m'|m) (™ + l|"i"> + E*s (m'|m + 1) (m|m")) . 

Because states with rh ^ ifi' are orthogonal, the sum 
reduces to a compact expression 
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FIG. 9: (a) Matrix elements of charge operator A'' (in units of Cooper pair charge 2e) for transitions g -k^ e, g -k^ f , and e /. 
(b) Matrix elements of phase operator (in radians) for same transitions. 
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FIG. 10: Top row: Effective potential for the "black-sheep" junction phase and a sketch of the three lowest levels for three 
values of mext. Bottom rows: Overlap between 27r-shifted fiuxonium wavefunctions in states g, e, and /, for the three values of 
nicxt. Note that for states g and e the overlap is clearly maximal when |mext| = 1/2. 



where the overlap function Fap{rnc:^t) is defined by sian di stribution with zero meanand standard deviation 

Eq. The flux dependence of v^^^ comes from the a = \J {Jie[Es\Y, 

flux-dependcncc of the qubit state wavefunction, several 
examples are illustrated in Fig. [lO]and in Fig. [TTJ 



In order to convert the shift v\ 



, Eq. (|B2|), into the 
linewidth Si^ap, Eq. let us recall that 



AT 

E 



Es- exp {i2iiQj/2e) 



(B3) 



with Qj being random variables with a spread of val- 
ues comparable to e, and the sum running over all array 
junctions. According to the central-limit theorem, in the 
limit of large N , the quantity Re[i?s] obeys to the gaus- 



P{0 < Re[£;5] < x} = —L^ exp(-a;V2cr2) . (B4) 



Assuming the array junctions to be approximately iden- 



tical, Eq 



EsA, we get {Re[Es] 



1 p2 



and then readily compute the linewidth Svap (defined 



(^i^'^) ) inhomogeneous broadening to be 

given by Eq. (|6]). If charges Qj vary slowly compared 
to the duration of a single Ramsey fringe experiment 
(of order 10 /xs), the decaying envelope of the of Ram- 
sey fringes is given by the absolute value of the charac- 
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Flux bias m^^, = <t>^J<S>^ 



FIG. 11: Flux-dependent part of the inverse linewidth due 
to coherent quantum phase-slips and the Aharonov-Casher 
effect in the array. Black line corresponds to the model A 
prediction for the g -f-*- e transition, purple and yellow lines 
present results of model B for the g o e and g f transi- 
tions, respectively. Note that close to |mcxt| = 1/2 the two 
models provide the same result for the lowest transition. 



teristic function of the distribution (B4|. We therefore 



with the flux-dependent dephas- 



ing time jj^Q^^ of the a — > /3 transition due to the co- 
herent quantum phase-slips given by 



me.t= 0.000 



m, 



ext 



0.005 



FIG. 12: Sketch of the lowest three energy levels of fluxonium 
qubit and structure of the corresponding relaxation rates in 
close proximity to rriext = 0. Vanishing of the phase matrix 
element ip^^ at rriext = combined with the particular fre- 



quency dependence of -Rintrinsic 
T/^ at zero flux bias. 



iuS) explains the sweet-spot in 



find that Ramsey fringe envelope is given by a gaussian rule' 



1/Tf^ 



47r- 



Rq 



, coth ■ 



2kBT 



Here, 



(B7) 



(B5) 



We used the expression (B5) to produce theory plots in 
Fig.i 



2. Energy relaxation of fluxonium transitions (Ti 
processes) 

Transitions e — > g (g — > e). Energy relaxation of the 
qubit, taking into account the finite temperature of the 
sample, an important effect around moxt = 1/2, takes 
place at the rate Fi — Tg^^ + ^e^g, where Ta-^|3 is the 
rate at which the qubit makes a transition from states a 
to /3. The relaxation time is defined as T^^ = l/Fi and is 
given in terms of the black-sheep junction phase matrix 
element ip^g (Fig. pb) and effective frequency-dependent 
parallel resistance 7toff(w) shunting the black-sheep junc- 
tion. We may decompose this resistance into the compo- 
nents i?intrinsic('^) and i?pui.coii (i^) : Originating from the 
spontaneous emission into the dissipative bath associated 
with the qubit circuit, and into the measurement appa- 
ratus, respectively: 



7->-i _ n-l 

-^cff ^ -'^intrinsic 



(B6) 



The latter dissipation mechanism comes from the 
resonator-filtered 50 fl environment of the measurement 
circuit (Fig. [7^-b) and is called Purcell effect . The re- 
sulting formula for T^^ is obtained from Fermi's golden 



R 



Pureed 



11 (^) = 



— 1 i^^ociCc 



(cj/cjo) 



4cot2(|w/wo) 



(B8) 



is the real part of the admittance of the electrical circuit 
connected to the black-sheep junction. In Fig. [6^ this is 
a circuit connected to the two terminals of the element 
Zl"^(w) with a pair of coupling capacitances 2Cc (the 
capacitances 2Cc are also shown as the interdigitated ca- 
pacitances in Fig. [l]i and Fig. The above expres- 
sions work as long as |1 — iz/i^ol ^ Qcxt- Note that 
^Pmccii ^ 1/C'c, and that the value of Cc in our circuit 
is of order 1 fF, making the Purcell contribution smaller 
than in transmon qubits (Cc ~ 20 fF) by more than two 
orders of magnitude. For the parameters of our sample, 
the Purcell contribution becomes negligible as soon as 
I'D ~ i^eg > 300 MHz. Therefore, energy relaxation in 
our qubit is mostly intrinsic. Using the reasoning in the 
main text, we express i?intriiisic through the effective (and 
also frequency-dependent) shunting resistances Rj of the 
junctions (Fig. [TJj): 



R- 



intrinsic 



— Rn 



N 



(B9) 



Since the area of the black-sheep junction is only a factor 
of 6 — 7 smaller than that of the array junctions, and 
assuming that RJ^ is proportional to the area, it is likely 
that for large N, only the black-sheep junction (j — 0) 

contributes, i?intrinsic ~ Ro- 
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Relaxation f ^ g. In the vicinity of rrioxt = 0, we 
have Te^g ^ r/->e, In the two-step transition 

f e g the step, / — e, is the bottleneck (Fig. |9]d, 
12 1. Therefore, the relaxation time associated with the 
f ^ g transition can be written as 1/t/^ ~ T f^e+^ f^g, 
and, neglecting the Purcell contribution, we get 



An- 



Rq 



+ 477 



Rq I |2 



(BIO) 



R 



We dropped the temperature-dependent factor in the 
first term because, in our experiment, in the vicinity of 
Wext = 0, the transition energy are such that hvfe ~ 
fcsT ^ hv fg. Once ksT ~ hv fg^ the relaxation process 
becomes more complicated. 



3. Common dephasing mechanisms 



Noise in mext- We establish a higher bound on its 
amplitude, assuming the noise is "l/f {Sm^^t — ^^txtl^^ 
where v is the frequency in Hz), from pure dephasing 
time at moxt = 0.2, where the T2 is maximal (Fig. 3a) 
but Aharonov-Casher contribution is suppressed. For a 
"l/f flux noise the pure dephasing time (gaussian decay 
of echo signal) is given bjEll 



^moxt X dl'eg("^cxt)/dmext 



(Bll) 



We estimate the dephasing time by subtracting the 
decay constant \/T2 of the echo signal (which was 
nearly exponential), from the separately measured l/2ri: 
ymext ^ _ i/2ri)-i > 35 /iS. Given that 

dj^eg(wext)/dTOext — 19.2 GHz (Fig. [2jD-c) we thus ex- 
tract the "l/F flux noise amplitude to be (Jwoxt < 
2 X 10-6. 



Noise in Ej = Ej^ . According to a previous stud; 
this noise is believed to be "l/f with the spectral density 



Sej — SE'j/v, where 5Ej (x Ej. To first order, the de- 
phasing time T2 due to fluctuating Ej is proportional to 
(/idz^eg(moxt)/d-Ej)~^. Remarkably, present fluxonium 
circuit is supposed to be insensitive to the noise in Ej 
for some specific value of mext — 0.43 (Fig. 13 1). Over- 



all theoretical non-monotonic dependence of T2 on flux 
makes the noise in Ej completely incompatible with the 
data. For instance, if the Ej noise limits dephasing time 
at mext — 1/2 to the measured value of 250 ns, then it 
should also limit it to a similar number at mext — 1/4, 
where we measure the dephasing time one order of mag- 
nitude longer. By analogy with the flux noise relation 
(Bll ), we can estimate 5Ej from the similar relation 



1/Tf ^ 6Ej X dz.eg(mext)/d-B, 



(B12) 



Assuming that at mext = 0.2 the pure dephasing is en- 



tirely due to "l/f Ej noise, we substitute T^-' > 35 /iS 
and hdveg{rnext)/<iEj 



0.1 (Fig. 13 5,), we extract a con- 
-^1 



servative estimate 6Ej < 3 x 10 

Noise in E^. Fluctuations in the Josephson ener- 
gies Ej^^g of the array junctions would result in a noisy 

inductance (noisy El) such that Se^ — ^-E-l/^- ^1" 
ready from the considerations of model A it is clear that 
d;/ge(mext)/di?L is minimal at mext = 1/2 and maxi- 
mal at mext = 0, see Eq. PI and Fig. [If. Calculating 
this derivative numerically (Fig. 13 d) we conclude that 



the theoretical flux-dependence of the E^ noise indeed 
is completely inconsistent with our dephasing data. We 
extract the higher bound on the E^ noise, assuming it is 
1 /f , from the pure dephasing time of the g ^ f transition 
measured at mext = 0, where flux, Ej, and Aharonov- 
Casher dephasing effects are minimal (Fig. [2j:, Fig. 1 
and Fig.fisk). Using, as usual, l/T^^ = I/T2 - l/2r/ 
the relation 



l/T^^ = SEl X diyfg{m,^t)/dEL 



(B13) 



the observed values > 50 /^s, and hdv f g{me^t) / dE ^ 
16 (Fig.fTsl)) at mext = 0, we find 5El < 3 x IO^^^l- 



16 




Elapsed time (Ks) 

FIG. 14: Measurement of T\ (blue trace), T2 (red trace) and T2 (green trace) as a function of time for the g e transition at 
jwicxt] =0.2. A set of the three data points was acquired continuously every 8 seconds. 



4. Time-stability of the decoherence times 



Fig. 14 shows that the measured values of Ti may show 
fluctuations by as much as 50% over time. The data is 
taken for the lowest g -f-)- e transition for rricxt = 0-2 
and is typical for any other bias. The origin of these 
Ti fluctuations is unknown. Fortunately, the dephasing 



times T2 , measured in a Ramsey fringe experiment, being 
significantly lower than 2Ti, are almost unaffected by 
the fluctuations in Ti. We also note that the times T2, 
measured in an echo experiment, follow the fluctuations 
of Ti , confirming that the value of T2 is indeed close to 
2Ti. 
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Readout resonator strips width (measured) 


15 /im 


Readout resonator strips separation (measured) 


4 /im 


Readout resonator wave impedance Zoo (inferred) 


80 n 


Readout resonator resonance frequency uq (fit) 


8.175 GHz 


Readout resonator external quality factor Qoxt (measured) 


400 


Readout resonator internal quality factor Qint (measured) 


2000 - 5000 


Black-sheep junction dimensions (nominal) 


0.35 /imxO.2 nm 


Array junction dimensions (nominal) 


2 pmxO.2 /im 


Array inductive energy El ~ {h/2e)^ /L (fit) 


0.525 GHz 


Array inductance L (inferred) 


300 nH 


Black-sheep junction Josephson energy Ej^^g = Ej (fit) 


8.9 GHz 


Black-sheep junction Coulomb energy Ec^^g = Ec = e^/2Cj (fit) 


2.5 GHz 


Number of array junctions A'^ (nominal) 


43 


Array junction Josephson energy Ej^tQ = NEl (inferred) 


22.5 GHz 


Array junction Coulomb energy Ec^^^ (inferred) 


0.85 - 1 GHz 


Qubit-eavity coupling constant g (fit) 


181 MHz 


Qubit-eavity coupling capacitance Cc (inferred from g and Zoo) 


0.8 fF 


Black-sheep junction phase-slip energy \Esb\ (measured, inferred) 


369 MHz 


Array junction phase-slip energy Esa (inferred from L, Ej^^^, Ec^^q) 


50 - 250 kHz 


Array junction phase-slip energy Esa (inferred from T2*(mext = 1/2)) 


130 kHz 



TABLE I: Device parameters 



